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ABSTRACT 


An  N-component  parallel  system  is  subjected  to  a  known  load  program. 
As  time  passes,  components  fail  in  a  random  manner  which  depends  on 
their  individual  load  histories.  At  any  time  t  ,  the  surviving  com¬ 
ponents  share  the  total  load  according  to  some  rule.  The  system's 
lifetime  distribution  is  studied  under  various  breakdown  rules. 

Under  the  linear  breakdown  rule  it  is  shown  that  if  the  load  program 
is  increasing  the  system  lifetime  is  IFR.  Using  the  notion  of  Schur 
convexity,  stochastic  comparison  of  different  systems  is  obtained. 

It  is  also  shown  that  the  system  failure  time  is  asymptotically  nor¬ 
mally  distributed  as  the  number  of  components  grows  large.  All  these 
results  hold  under  various  load  sharing  rules,  in  fact,  one  can  prove 
that  the  system  lifetime  distribution  is  invariant  under  different 
load  sharing  rules. 

For  a  more  general  breakdown  rule  only  the  equal  load  sharing  is 
considered.  The  asymptotic  distribution  of  the  system  lifetime  is 
shown  to  be  normal.  The  asymptotic  mean  and  variance  are  derived. 
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CHAPTER  1 


INTRODUCTION  AND  SUMMARY 

This  report  is  concerned  with  the  study  of  parallel  systems  whose 
components  are  subjected  to  a  certain  nonnegative  time  dependent  load. 

Consider  a  parallel  system  of  n  components.  The  system  sustains 
a  certain  load  (damage,  wear)  which  varies  in  time,  and  thus  each 
functioning  component  is  exposed  to  a  certain  fraction  of  this  load. 

As  time  passes,  components  fail  in  a  random  manner  which  depends  on 
the  amount  of  damage  they  have  been  exposed  to.  As  the  single  compo¬ 
nents  fail  one  by  one,  each  of  the  surviving  components  is  exposed  to 
an  increasing  share  of  the  total  load.  This  model  is  known  in  the 
literature  as  the  "load  sharing  model"  and  has  been  treated  previously 
by  Coleman  (1958) ,  Birnbaum-Saunders  (1958)  and  more  recently  by 
Phoenix  (1978).  Coleman  adapted  this  model  in  connection  with  the 
strength  behavior  of  fiber  bundles  and  his  theoretical  results  were 
remarkably  consistent  with  the  experimental  behavior  of  a  variety  of 
structural  materials.  Birnbaum  and  Saunders  used  the  model  to  study 
the  "Fatigue"  of  materials  under  conditions  involving  dynamic  loading. 
They  also  reported  a  very  good  agreement  with  empirical  data.  Phoenix 
(1978)  adapted  Coleman's  assumptions  about  the  stochastic  behavior  of  the 
single  component  under  a  given  load  program.  He  made  use  of  Stig]er 
(1974)  paper  to  show  that  the  system  failure  time  is  asymptotically 
normally  distributed  as  the  number  of  components  grows  large  and  thus 
generalized  Coleman's  results  which  dealt  with  the  calculation  of  the 


asymptotic  mean  time  to  failure.  All  previous  study  assumes  that  all 
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Che  components  are  idenCical,  that  is  each  individual  component  has 
the  same  stochastic  behavior  under  a  same  load  program.  Furthermore, 
the  surviving  components  at  any  time  t  are  assumed  to  share  the  total 
load  equally. 

In  this  report  we  generalize  these  assumptions.  We  allow  the 
components  to  be  nonidentical  and  for  some  special  case  to  share  the 
total  load  in  a  fairly  general  way. 

In  Chapter  2  we  adapt  the  linear  breakdown  rule  that  was  essentially 
used  by  Birnbaum-Saunders.  Loosely  speaking,  if  a  single  component  is 
subjected  to  a  time  dependent  load  program,  the  failure  rate  of  the  com¬ 
ponent  at  any  time  t  is  linear  in  the  load  at  t  .  Birnbaum-Saunders 
were  interested  in  certain  statistical  aspects  of  the  system  (e.g., 
estimating  the  number  of  components  in  the  system  and  the  deterioration 
factor  of  the  individual  component).  We  focus  on  the  stochastic 
characteristics  of  the  system.  Among  other  things  we  prove  that  the 
components  lifetimes  are  stochastically  associated,  that  if  the  system's 
load  program  is  nondecreasing,  its  failure  time  is  IFR.  Using  the 
notion  of  Schur  convexity,  we  obtain  stochastic  bounds  on  the  system 
lifetime  distribution.  Asymptotic  results  are  also  available  as  the 
number  of  components  in  the  system  grows  large.  All  these  results  are 
shown  to  hold  under  a  fairly  general  load  sharing  rule.  In  fact,  we 
prove  that  the  distribution  of  time  to  system's  failure  is  invariant 
under  the  load  sharing  rule. 

In  Chapter  3  we  adapt  Phoenix  model  and  generalize  it  to  the 
nonidentical  component  case.  We  obtain  asymptotic  results  similar  to 
those  in  Phoenix  (1978)  and  also  point  out  an  unjustified  statement  in 


his  paper  and  give  our  version  for  it. 
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CHAPTER  2 

THE  LINEAR  LAW  BREAKDOWN  RULE 

Consider  a  system  of  N  components  in  parallel  which  are  being 
subjected  to  a  nonnegative  time  dependent  load  Ls(*)  •  At  any  time 
t  ,  all  the  surviving  components  share  the  load  according  to  some  rule. 
In  the  first  three  sections  of  this  chapter,  we  assume  that  they  share 
the  load  equally.  This  assumption  is  removed  in  Section  2.4,  where 
we  consider  a  more  general  rule.  Thus  for  the  time  being,  if  i 
(i  <  N)  components  have  failed  prior  to  time  t  ,  the  actual  load  on 
each  of  the  N  -  i  surviving  components  at  t  is  Lg(t)/N  -  i  . 

Component  i  (i  =  1,2,  . . . ,  N)  subjected  to  a  known  load 
program  1(«)  has  a  random  failure  time  T  .  Throughout  this 
chapter  we  assume  that  under  l(-)  the  distribution  of  T  is  of 
the  form: 

t 

-Bi  /  l(u)du 
F1(t  |  1)  =  1  -  e  0 

where  is  a  strictly  positive  constant  which  characterizes  the 

durability  of  component  i  . 

If  we  define 

fi(t  I  1}  =  dt  FiU  I  1}  and 

fi(t  |  1) 

*  1  -  F1(t  |  1)  C  -  ° 

then  it  is  easy  to  check  that  A^(t  |  1)  -  J  1 (t)  . 


^(t  |  1) 
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A^(t  |  1)  is  the  failure  rate  function  of  component  i  under  the  known 
load  program  1(*)  .  It  can  be  interpreted  as  the  instantaneous  prob¬ 
ability  of  a  t  years  old  component  i  to  fail,  that  is,  component 
i  having  not  failed  prior  to  time  t  and  sustaining  load  l(t)  will 
fail  during  the  interval  (t  ,  t  +  At)  with  probability  Sil(t)At  +  o(At) 
where  o(At)  denotes  remainder  terms  of  order  less  than  At  as  At  -*■  0  . 
For  the  obvious  reason  we  say  that  the  component  obeys  the  linear  law 
breakdown  rule.  As  a  system,  the  components  lifetimes  are,  therefore, 
strongly  dependent  but  only  through  their  failure  rate  functions.  For 
instance,  if  at  t  there  are  j  (1  <  j  <_  N)  functioning  components  and 
a  component  fails,  then  under  the  equal  load  sharing  rule  the  load  that 
each  of  the  surviving  components  sustain  changes  instantaneously  from 

L  (t  )/j  at  t  to  L  (t)/j  -  1  at  t  . 
s  s 

2.1  The  Joint  Distribution  of  Components  Lifetimes 

We  denote  by  T^,T2>  •••,  T^  the  lifetimes  of  components  1,2, 

...,  N  respectively  and  by  <_  T^)  —  ’  •  •  £  T(n)  their  order 

statistics.  We  assume  that  the  system  is  being  subjected  to  a  non¬ 
negative  load  program  Ls(t)  and  that  each  component  obeys  linear 
law  breakdown  rule.  Each  of  the  components  is  being  characterized 
by  its  proportional  factor  8  (Birnbaum,  Saunders  (1958)  called  it 
the  deterioration  factor)  which  is  assumed  to  be  strictly  positive 
number.  The  system  itself  is  being  characterized  by  its  load  program 

L  (•)  and  by  the  vector 
s 

£  -  (81,82,  ....  8n)  • 

The  above  assumptions  lead  to  the  following  facts: 


5 


(2.1)  For  any  i  (i  =  1,2,  . . . ,  N)  , 

L  (t) 

P(Ti  e  (t  ,  t  +  At]  I  T(1)  >  t)  =  ei  — jj —  At  +  o  (At) 
as  At  -*•  0  . 

(2.2)  If  0  <_  t^  <_  t^  £  . .  •  £  t^  <_  t  ,  k  <  j  <_  N  ,  and  {n^.n^ ,  •  •  •  ,  n^} 
is  an  arbitrary  permutation  of  {1,2,  ...,  N}  ,  then: 

P(V.  e  '  +  “1  I  \  '  V  •••■  \  ■  ck  •  T(k+1)  *  c)  ' 

L  (t) 

6n.  At  +  0(At> 

J 

as  At  -*■  0  . 

Facts  (2.1),  (2.2)  tell  us  that  if  exactly  k  (k  =  0,  . . . ,  N  -  1) 
components  failed  by  time  t  and  component  j  is  not  among  them, 
then  the  failure  rate  of  j  at  t  is 

L  (t) 

6j  N  -  k  * 

From  (2.2)  we  also  learn  that  the  order  in  which  the  components  failed 
does  not  affect  the  failure  rate  of  component  j  ,  as  long  as  it  is 
given  that  n  -  j  survived  t  and  j  is  among  them. 

Lemma  2.1: 


P (T(1)  >  t)  =  e 


e  +  ...  +  e  t 

—  /  L  (u)du 
0  s 


N 
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Proof : 


From  (2.1),  it  is  obvious  that  the  failure  rate  function  of  is 


Ls(t)  N 


-  ■ 


Q.E.D 


Lemma  2.2: 

Let  0  <  t1  <_  t2  £  . . .  £  tk  <_  t  ,  k  <  N  and  let  {n^n,,  ,  .  .  .  ,  r^} 
be  an  arbitrary  permutation  of  {1,2,  ...,  N}  .  Then  for  x  >  0  we  have 


P\T(k+l)  +  "  Ck  *  V+D  =• 


6  +  . . .  +  6  t+x 


exp 


(  p 

)  k+1  °N 


N  -  k 


I 


L  (u)du;  . 

s  i 


Proof: 


Follows  from  (2.2)  in  a  straightforward  way. 


Lemma  2.3: 


Q.E.D 


(t^}  ,  {n^}  as  in  Lemma  2.2  and  let  d_t  denote  the  infinitesimal 
interval  (t  ,  t  +  dt)  ,  then  for  x  >  0  : 


'{T  G  dt  , 
\  \  ~ 


(k+1) 


>  t  +  x  I  T  =  t. . T  =  t,.  ,  ,  T 


°1 


nk_x  k-1  •  *(k) 


)  = 


N""'k  +  1  Ls(t)  exp 


I  6" 

r 


k+1 


"n 


N  -  k 


t+x  I 

J*  Lg(u)du 


dt 
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Proof: 


-It 
\  nv 


p|\e^  ’  T(k+D  >t  +  *\\mtV  •  —  TVl  =  Ck-1  •  T(k)  >  C1  - 


PIT,,  >  t  +  x  T  =  tlt  T  =  t.  ,  ,  T  =  t.  .  ,  T  =  t  , 

\  (k+1)  1  n1  1’  nk_1  k-1  nk  k-1  nR 


T  >  t 

(k+l) 


}  •  PIT  €  dt  I  T  =  t,,  ....  T  =  t.  .  ,  T  >  t)  . 

)  \  nk  ~  nl  1  Vi  k_1  (k)  / 


Combine  (2.2)  and  Lemma  2.2  and  the  result  follows. 


We  are  now  ready  to  state  the  following  theorem: 


Theorem  2.1: 


Q.E.D. 


Consider  the  simplex  0  <  <  . . .  <  and  let  {n^.n^,  . 

a  permutation  of  {1,2,  ...,  N}  ,  then 


>(T  G  dt.  ,T 

\ni  — 1  n 


G  dt^2>  •••>  Tn  e 

2  N 


...  „N)  be 


N 

n  0.L  (t.)dt. 

>i  i  s  j  j 

N! 


•  exp  - 


8  +  .  . .  +  6  1 


“N 


J 


L  (u)du 
s 


8  +  ...  +  8  2 


^2 _ y  f 

N  -  1  J 


L  (u) du  +  ...  +8 
S 

C1  tN-l 


N 

/ 


Lg (u)du 
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Proof : 


P  T  e  dt.  ,  . 
\  nl  — 1 

01 

H 

P|T  G  dt,,  1 

t,  e  dt, , 

.  . .  ,  T 

K  -N  1 

i 

l  — i’ 

Vi 

x  PlT  €  dt„ 

,  ,  T  > 

tH  1  t  € 

\  vi  -*■ 

-i 

N  1  n^ 

T(N-1)  >  CN-1 

) 

-  p(\€^i  ■ 

T(2)  >  C2 

1  T(i)  >  c: 

c») 


‘i-  V2  ' 


«  P(T(1)  »  tx)  ■ 


Combine  Lemmas  2.1,  2.2,  2.3  and  the  result  follows. 


Corollary  2.1: 

If  8^  =  =  •  •  •  =  ^  ~  ®  ’  then  ttle  density  function  of 

given  by 


£-  (t)  •  SNL  (t) 

T<») 


,  N-l 


Lg(u)du  J 


-6  /  L  (u)du 

J  Q 


(N-l)! 


Proof : 

By  Theorem  2.1,  we  have  that  on  0  <  t,  <  . . .  <  t„ 

IN 


P  (T 


(1) 


e  dt_1  > 


(N) 


£  dt^) 


N 

n  L  (t.)Bdt. 
j=l  S  ^  J 


-6 

e 


t 

/  L  (u)du 
0  s 


Q.E.D. 


is 
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Integration  yields  the  result. 


Q.E.D. 


In  fact,  this  result  can  be  obtained  directly  from  the  simple 
observation  that  if  the  0's  are  equal  then 


P(T(i)  E~  1  T(l)’  T(i-1)  ’  T(i)  >  C) 


SLs(t)dt  on  T(._1)  <  t 


on  T(._1)  >  t  . 


Thus  the  ordered  sequence  of  failure  times  T.,.  <  T.„x  <  ...  <  T,,,. 

(1)  —  (.2)  —  —  (N) 

behaves  like  a  sequence  of  "events"  taken  from  a  nonhomogeneous  Poisson 
process  with  intensity  function  0Ls(*)  ,  and 


P(T(N)  <  =  Pr  {there  are  at  least  N  events  of  the  Poisson 

in  (0 , t  ]  } 


=  1  - 


N  /  t  V- j 

l  8  /  L  (u)  /j! 

j=0  '  0  s  > 


-0  /  L  (u)du 
0  s 


and  —  P(T^^  t)  yields  the  result. 

Let  us  now  define  the  random  process  L ( t ) 


\(t)/N  0  <  t  <  T(1) 


(2.3)  L( t) 


Ls(t) 
iN  -  i 


T(i)  -  C  <  T(i+1)  1 


=  1,2 . N  -  1 


Ls(t) 


t  >  T 


(N) 


r 
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1 


and  let 

T. 

i 

(2.4)  Qi  =  J  L( u)du  ,  i  =  1,2 . N  . 

0 

Thus  intuitively,  is  the  total  load  sustained  by  component  i 

during  its  lifetime. 

Theorem  2.2: 

The  random  variables  Q.,Q0,  •••,  Q,,  ate  independent  exponentially 

X  Z  N 

distributed  with  rates  •••»  respectively. 

Proof: 

The  joint  density  function  of  the  Q’s  can  be  written  in  terms  of 
the  joint  density  of  the  T's  : 

•••»  =  •••>  tjj)||j|| 

where  t±  =  ti(q1,  . . . ,  qN)  ,  i  =  1,2,  . . . ,  N  and  |  | J|  j  = 

is  the  Jacobian  of  ^  with  respect  to  j:  . 

Now  consider  the  simples  0  <  q^  <  q£  <  ...  <  q^  and  let  t^,t2>  •  •., 

t.T  be  such  that 
N 


3(q1.  • 

t*N> 

a(tr  . 

. • ,  tN) 
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L  (u) 


r  ^ 

h  *  J  — 


du 


L  (u) 


q2  =  q 


/• 

1  J  N  -  1 


du 


’  Vl  +  / 


L  (u)du  . 
s 


N-l 


Clearly  <  ...  <  and  the  joint  density  function  of  Q^, 

on  the  above  simplex  is 


f£(ql’  •*'»  V 


^  6iLs(ti(q1,  ....  qN))e 
_ 


■6iqi 


i-l 


By  simple  calculus  we  also  have 


/ 


Ls<Cl> 


w 

N-l 


0 


.L  (tN) 
s  N  i 


and  hence 


N 

n  l  (t  ) 
i=l  S 
N! 


Thus  on 


0  <  <  q2  <  •••  <  %  *  c^e  density  on  the  Q's  is 


qN}  = 


N 

n 

j=i 


8 .  e 
1 


-8.q. 
J  J 


It  is  fairly  easy  to  check  this  argument  holds  for  an  arbitrary  simplex 
of  the  form  0  <  <  ...  <  q^  ,  where  {n^,n2»  ....  n^}  is  a 

permutation  of  {1,2,  ...»  N}  .  Thus  the  joint  density  function  of 

Q1,Q2.  •••>  Qn  is  given  by 


^  ~®i^i 

fn  n  (q. ,  •••,  qN)  *  n  B.e  ,  q  >  0  ,  i  »  1,2,  ...,  N  . 

^1 ’  *  *  * ’  1  i*l 

Q.E.D. 

The  importance  of  this  theorem  is  that  it  establishes  the  existence  of 

a  transformation  *  ^1’  ^1’  Qjj)  wbich  under  a 

s  s 

given  load  program  bs(”)  maps  the  random  vector  with  dependent  com¬ 
ponents  into  with  independent  components.  This  fact  enables  us  to 
analyze  most  of  the  parameters  of  the  system  in  terms  of  independent  random 
variables.  Indeed,  throughout  the  next  two  sections  we  use  this  fact  to 
find  conditions  under  which  the  system  lifetime  is  IFR,  obtain  bounds 


for  the  system's  reliability  function  as  well  as  to  establish  the 
asymptotic  distribution  of  the  system's  lifetime  as  N  ,  the  number 
of  components,  increases  to  infinity. 

We  start  with  proving  the  rather  intuitive  fact  that  the  life¬ 
times  of  the  individual  components  (T^,  T^)  are  associated,  but 

first  some  necessary  preliminaries. 

Definition: 

We  say  that  the  random  variables  Z,  ,Z. . Z  are  associated  if 

i  l  n 

Cov  ( f (Z) , g ( Z) )  j>  0  for  all  nondecreasing  functions  (f,g)  for  which 

the  covariance  is  well  defined.  We  use  Z  for  (Z,  ,  ...,  Z  )  . 

In 

For  completeness  we  start  with  the  following  properties  of  associa¬ 
tion: 

(P^)  Any  subset  of  associated  random  variables  are  associated. 

(P2)  If  two  sets  of  associated  random  variables  are  independent 
of  one  another,  then  their  union  is  a  set  of  associated 
random  variables. 

(P^)  The  set  consisting  of  a  single  random  variable  is  associated. 
(P^)  Nondecreasing  functions  of  associated  random  variables  are 
associated. 

(P^)  Independent  random  variables  are  associated. 

From  the  way  :  (T^,  ...,  T^)  -*■  (Q^,  . ..,  Q  )  is  constructed, 

s  L 

it  is  clear  that  if  the  load  program  Lg(*)  is  strictly  positive 

for  each  t  ,  ^  is  one  to  one  and  hence  is  also  well 

s  s 

Now  in  order  to  prove  that  (T^^,  ...,  TN)  are  indeed 


defined. 
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associated,  it  suffices  to  show  that  is  nondecreasing,  that  is, 

s 

T^(Q^»  •••»  Qjq )  >^2  ^i*  ■■■>  ’  * "  *  *  *  •  •  • »  )  are  all  non¬ 

decreasing. 

Lemma  2.4: 

If  L  (•)  >  0  ,  then  t|i  ^  is  well  defined  and  increasing 
S  Li 

s 

functional. 


Proof: 

The  first  part  follows  easily  from  the  way  is  constructed. 

s 

To  show  that  :  (q  ,  ....  q^)  -*■  (t^,  ...,  tN>  is  increasing  func- 

s 

tional,  denote  by  q^  £q(2)  £  ...  iq(N)  and  C(X)  1  t(2)  C(N) 

the  increasing  rearrangement  of  q^,  . . . ,  q^  ,  t^ ,  . . . ,  t^  respectively. 

Under  •••»  ^  (2)  ^1  ’  ' ' '  ’  ’  ^(N)^l’  ‘ 

are  such  that 


l(D 


f  Lidl 

J  N 

0 


du 


‘(2) 


’(I) 

J 


Ls(u) 


du  + 


T  ^ 

J  N  -  1 


du 


(1) 


C(l)  T  (U) 

■  J  -V 

o 


+ 


t 


N 


J 


N-l 


L  (u)du  . 
s 
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Now  fix  all  q's  and  increase  q..,  by  the  infinitesimal  increment 

(i) 

dq^^  .  It  is  easy  to  check  that  for  j  <  i  ,  ^  remains  unchanged. 

For  t^  we  have 


<“>(1)  ■  d‘(i>  ■ 


To  find  dt...,.  note  that 
(i+l) 


Ls(t(i+1)) 

0  ~  dq(i+l)  "  dq(i)  N  -  i  dC (i+l) 


Ls(t(i)) 
N  -  i 


dt 


(i) 


or 


dt 


(i+l) 


Ls(t(i+l))  N 


dt 


(i) 


i+l 


and  indeed  for  k  >  1 


dt 


(i+k) 


t-s<t(t))  dt 
Ls(t(i+k))  N  *  1  +  1 


iil 


1,2, 


which  completes  the  proof. 


Q. E.D. 


Corollary  2.2: 

The  components  lifetimes  T^,  . . . ,  T^  are  associated. 

Next  we  show  that  under  increasing  load  program  the  system  lifetime 
is  IFR. 
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Definition: 

A  nonnegative  random  variable  X  with  distribution  F(*)  is  said 

to  be  IFR  (or  equivalently  is  said  to  have  an  increasing  failure  rate 

distribution)  if  log  (1  -  F(x))  is  concave  where  finite. 

If  F  also  possess  a  density  f(-)  ,  then  it  is  easy  to  check  that 

the  above  is  equivalent  to  saying  that  X(t)  =  7 — ^  ^  is  non- 

1  -  r  (,  t ; 

decreasing  on  {t  :  F(t)  <  1}  . 

Lemma  2.5: 

Let  g  be  an  increasing  differentiable  function  and  let  (X,Y) 
be  two  random  variables  such  that  Y  *  g(X)  and  suppose  Y  has  a 
failure  rate  function  \  ,  then  the  failure  rate  function  of  X 

is  given  by  \j.(t)  »  A^.CgCt)  )g'  (t)  . 


Proof : 


f  (t)  f  (g(t))g’(t) 

1  -  Fx(t)  “  1  -  FY(g(t))  =  \(g(t))g'(t)  . 


Q.E.D. 


Theorem  2.3: 

If  the  system  load  program  Lg(*)  is  nondecreasing,  then  the  system 
lifetime  is  IFR. 

Proof : 


Let  Ql,Q2, 


Q„  be  as  defined  in  (2.4),  then  check  that 

1*1 
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N 

Now  since  the  Q's  are  independent  IFR,  so  is  £  Q  and  by  Lemma  2.5 

i-l  1 


we  have 


(N) 


(t) 


XEqH  Ls(u)|Ls(C> 


Hence  X_  is  nondecreasing. 
(N) 


Q.E.D. 


2.2  Schur  Functions  and  System  Reliability 

The  method  and  concept  of  majorization  and  Schur  functions  are  used 
in  this  section  to  make  stochastic  comparison  between  systems  with 
different  components  (that  is,  components  with  different  6-factor)  and  to 
develop  conservative  bounds  for  the  system  lifetimes.  For  the  sake  of 
completeness  we  review  some  of  the  definitions  and  tools  of  majorization 
and  Schur  functions: 

Given  a  vector  x  *  (x, ,  ....  x  )  let  xr,,  >  xr„,  >  ...  >  xr  . 

-  1  n  [1]  -  [2]  -  -  [n] 

denote  the  decreasing  rearrangement  of  x^,x9,  ...,  x  . 

Definition : 

A  vector  x  is  said  to  majorize  a  vector  x'  if 

j  j 

l  xm  1  l  x'm  J  =  1 . n  -  1 

i-l  1  J  i-l  1  J 


and 


n  n 

l  x[i]  l  x[i] 

i-l  1  J  i-l  1  J 


We  write  x  >  x' 
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Although  majorization  involves  comparison  of  vectors  of  order  n  , 
the  following  characterization  shows  that  we  need  only  consider  a  pair 
of  coordinates  at  a  time. 

Theorem  2.4:  (Hardy,  Littlewood,  Polya.  1952,  p.  47) 
m 

Let  x  x'  •  Then  x'  can  be  obtained  from  x  by  finite 
number  of  T- t rans format ions ,  where  a  T-transf ormation  changes  two 
coordinates  only,  and  in  the  following  way.  If  for  example  coordinates 
1  and  2  are  being  transformed,  then  T(y_)  *  (y^,y^,y^,  ....  y')  where 
max  { y , y 2 ^  >  max  {y^y^}  and  y^  +  y^  =  yj  +  y^  .  Furthermore  we 
m 

have  2.  1.  T(y)  . 

Majorization  represents  a  partial  ordering  in  lRn  .  A  Schur 
function  is  a  function  that  is  monotone  with  respect  to  this  partial 
ordering,  more  formally: 

Definition: 

m 

A  function  h(*)  satisfying  the  property  that  x.  >_  x'  implies 
h(x)  >_  (<)  h(x/)  is  called  Schur  convex  (Schur  concave)  function. 

A  useful  characterization  of  Schur  functions  is  provided  by  the 
fundamental  theorem  of  Ostrowski  (1964)  which  states: 

Theorem  2.5:  (Ostrowski) 

A  differentiable,  permutation  invariant  function  h  on  IRn  is 
Schur  convex  (Schur  concave)  iff 

(xi "  “  jt)  -  (-}  0  for  311  1  *  J  • 

'  1  j' 

We  use  this  to  prove  the  following: 
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Theorem  2.6: 

If  X^  ,  X2  are  independent  exponentially  distributed  random 

variables  with  rates  A^  ,  A^  respectively  and  Xj  ,  X^  are  independent 

exponentially  distributed  with  rates  Aj  ,  A^  respectively,  and  if 
m 

(A^.A^)  >_  ( A , A ^ )  1  then  for  any  t  >_  0  P(X^  +  X2  >_  t)  P(X^  +  t)  . 


Proof : 

By  using  the  Ostrowski  Theorem  it  is  not  difficult  show  that 


P(xL  +  x2>t) 


is  Schur  convex  in  (A, ,A_)  . 

J.  <- 


Q.E.D. 


Corollary  2.3: 

If  X^jX^,  X^  and  Y^,Y2,  -•*,  Y^  are  two  sets  of  independent 

exponential  random  variables  with  rates  (A,,  ....  A  )  and  (A* . A') 

1  n  1’  ’  n 

m 

respectively  and  if  A_  >_  A_'  ,  then  for  any  t  >_  0 

p<xi  +  *‘,1')  ip«i  +  'U'1  • 


Proof : 

Use  Theorems  2.4,  2.6  and  the  well  known  fact  that  if  ,  X9 

are  independent  and  X|  ,  X^  are  independent,  and  for  any  t  ,  PCX^  >_  t)  >. 
P(X’  >_  t)  ,  i  =  1,2,  then  for  any  t  :  P(X^  +  X2  >_  t)P(X^  +  X^  >_  t)  . 


Q.E.D. 
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Remark: 


Proschan  and  Sethuraman  (1974)  proved  Lemma  2.6  using  a  slightly 
different  approach. 


Corollary  2.4: 


Fix  t  >  0  ,  then 


PI  I  Lg(u)du  >  t 


is  Schur  convex  function  in  •••<  Pjj) 


Proof : 


Lg(u)du  -  Q:  +  ...  +  Qn 


which  are  independent  exponentially  distributed  with  rates  3^, 
respectively. 


.... 


Corollarv  2.5: 


Fix  t  >  0  ,  then  P(T.„.  >  t)  is  Schur  convex  function  in 

I.N; 

Pr32,  ....  Bn  • 


Proof : 


Follows  from  Corollary  2.4 


Q.E.D. 


Q .  E .  D  . 


The  following  theorem  provides  a  lower  bound  for  th*'  survival 


probability  of  a  system  under  an  arbitrary  load  program  L  (•)  • 

Theorem  2.7: 

Let 


j=l 


3. 

1 


then  for  any  t  >  0  , 


-3  /  L  (u)du 


P(T(N)  *  C)  -  6 


l  frl  l  (u)du 

j=o3‘U 


Proof : 

m 

First  note  that  (3-^  B^)  1  (3,3, 

follows  from  Corollaries  2.1,  2.5. 


3)  .  The  result 


Q.E. 


2.3  Asymptotic  Distribution  of  System  Lifetime 

In  this  section  we  obtain  the  asymptotic  distribution  of  the  system 
lifetime.  We  show  that  under  certain  conditions  the  system  lifetime 
normalized  properly  tends  in  distribution  of  the  standard  normal  as  the 
number  of  components  increases  to  infinity. 

Definition: 

Let  ...  be  a  sequence  of  independent  random  variables  such 
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that  the  distribution  of  is  and  E[X^]  =  0 

k  =  1,2,  .... 


Let 


2  2  2 

S  =  X.  +  .  . .  +  X  ,s  =  a,  +  ...  +  o 
n  1  n  n  1  n 

We  say  that  the  sequence  {X^}  satisfies  Lindeberg  Condition  iff  for 
each  e  >  0  : 


2 

s 

n 


n 


l 

k=l 


/ 

|x|>es 
1  *  n 


x“dF^(x)  +  0 


as  n 


00  # 


The  following  are  well  known  theorems  that  we  will  use  later  on: 


Theorem  2.8:  (Lindeberg-Feller) 

If  Lindeberg  Condition  holds,  the  distribution  of  the  normalized  sums 
S 

—  tends  to  the  standard  normal  as  n  °°  . 


n 


We  write 


N(0,1) 


for  convergence  in  distribution. 


n 


Theorem  2.9:  (Slutzky) 

If  X  tends  in  distribution  to  X  and  {A  }  ,  {B  }  are 
n  n  n 

two  sequences  of  random  variables  converging  in  probability  to  a  ,  b 

(constants)  respectively,  then  AX  +  B  tends  in  distribution  to 

n  n  n 


aX  +  b  . 
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f 


Corollary  2.6:  (Slutzky) 

Suppose  that  a^  is  a  sequence  of  constants  tending  to  00  ,  b 
is  fixed  number,  and  a  (X  -  b)  ^  X  .  Let  g  be  a  differentiable 
function,  then  an(g(Xn>  -  g(b))  +  g'(b)X  . 

Theorem  2.10: 

Let  X^jX^,  ...  be  independent  exponentially  distributed  random 

variables  with  mean  ...  respectively  and  suppose  that  there 

exist  0  <  L  <  U  constants  such  that  for  each  n  :  L  <  0  <  U  , 

—  n  — 

then  {X  -  0  }  satisfies  Lindeberg  Condition, 
n  n 

Proof : 

The  distribution  of  , 

-t/0 

Fk(t)  =  1  -  e  ,  t  >  0 


and  hence  for  A  >  0  , 

J  x2dFR(x)  =  e  '  k(A2  +  2A9k  +  202)  . 

X>A 

Now  since  0k  are  uniformly  bounded,  check  that  for  n  large  enough 


x's  ,  that  is, 


it  suffices  to  consider  only  the  right  tail  of  the 


-T  I  J  ~  0k)2dFk<x>  =J2  l  f  ~  V2dFk(x> 

sk*l  s  k=l  „ J 


a  x-0  I >es 

1  k 1  a 


n  x-0.  >es 
k  n 


\  l  j  (*  -  « 

s  k=l 

n  x-0,  >es 
k  n 

/  A,*w 


(x  -  0k)"dFk(x) 


n  x>es 


n  -es  /0 


-4-  £  e  n  k  (e2s2  +  2es  0 

s2  k-1  \  n  a 


,  +  20" 1 
k  k  j 


1  -c/n. 1/15 (  2  „2 

— j  ne  (e  nU 

nL  ' 


|e2nU2  +  2e/nU2  +  2U2j 


n  e-(eL/U)^(e2u2  +  2t\}2  /  ^  +  2U2/n) 


which  vanishes  as  n  -*•  00 


Corollary  2.7: 


Q.E.D. 


Suppose  that  there  exist  0  <  A  <  B  such  that  for  anv  n  ,  A  <  6  <  B  , 

—  n  — 


^  N 

J  Lg(u)du  -  l  1/Bj 
_o _ jsl 

(i 


N  (0 , 1)  as  N  ® 
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Corollary  2.8: 


Let  0  =  . 

k 


0  >  0 


Then 


Suppose  that  0,  =  0  +  o(—\  as  k 

k  U  / 


(N) 

J 


\ 


Lg (u) du 


-•/ 


0  2 
"  N(0 ,6 Z)  . 


Proof : 


:(«  \ 

I  Lg(u)du 


-  0  /  =  / N 


[  Lg(u)du 


hk  /I  \ 

k=i  l ,  s  I  k=i 


e/  . 


The  first  term  equals 


(N) 


N 


/n 


J  Lg  (u)  du  ~  l  \  /N  ^ 

i _ /  Z  e£/N’ 

N  I k=l  K 


(j. a>) 


X  U 


®  and 


J 
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As  N  <*> 


CL  4 


-*■  1  and  the  second  term  ^  ^  0, /N  -  6 


Lk=l 


and  by  Slutzky  theorem,  the  result  follows. 


Example  1; 


Suppose  that  L^t)  =  Lg  constant)  and  0^  =  0  + 


Q.E.D 


•tt) 


as  k  -*■ 


where  0  >  0  the  asymptotic  distribution  of  the  system  lifetime  T^  is 


^(“IF  •  r)  *  »(o.<9/Ls)2) 


as  N 


Example  2: 


Suppose  that  Lg(t)  =  6t  ,  6  >  0  ,  and  the  {0^}  are  as  in  Example  1 


Then 


&  U 


(!%-•) 


+  N(O,02)  . 


Use  now  Corollary  2.6  to  show  that: 


I(N)  ‘  YT  N  "  N(o»0/26)  as  N 


The  interesting  thing  about  this  example  is  that  for  linear  load  program 
the  asymptotic  mean  of  T^  ,  E[T^]  ~  whereas  the  asymptotic 
variance  is  constant. 


2/ 


A  close  study  of  Theorems  2.1  and  2.2  shows  that  the  equal  load 
sharing  rule  used  throughout  the  previous  sections  is  not  essential  for 
the  derivation  of  most  of  the  results  so  far.  It  seems  that  a  system  of 
components  which  obey  the  linear  law  breakdown  rule  is  invariant,  in 
some  sense,  with  respect  to  the  load  sharing  rule.  In  this  section  we 
assume  a  fairly  general  rule  and  prove  that  a  modified  version  of 
Theorems  2.1  and  2.2  holds.  Thus,  using  similar  arguments  as  before 
we  obtain  similar  results  (i.e.,  IFR,  Schur  convexity,  asymptotic  dis¬ 
tribution).  In  fact,  we  show  that  no  matter  how  the  total  load  is  being 
distributed  among  the  functioning  components,  the  lifetime  of  the  system 
remains  stochastically  the  same. 

Thus  for  instance,  in  analyzing  the  mechanical  breakdown  in  bundles 
of  fibers  the  "local  load  sharing"  rule  received  some  attention  in  the 
literature.  The  rule  here  is  that  every  time  a  fiber  fails,  the  adjacent 
fibers  take  on  its  load,  we  show  that  this  system  does  not  differ  much 
from  the  equal  load  sharing  system  or  from  a  system  that  allocates  the 
total  load  to  a  single  component  at  a  time. 


Again,  Dy  T^.T^,  . . . ,  T^  we  denote  the  lifetimes  of  components 
1,2,  ...,  N  respectively  and  T^  T^)  1  •••  T^  are  their  order 
statistics.  N(t)  is  the  number  of  components  that  failed  by  time 


t  (t  >_  0)  .  Let  Ht  denote  the  history  of  the  system  at  time  t  , 
that  is,  is  the  o-field  of  events  generated  by  the  events  of  the 


|Tn1ltl,T»2lt2’  -,\ltk  ’N(0  =kj 


where 


i  =  1,2,  ...,  k  ,  k  £  N  and  { n^ ,  . 


of  (1,2,  ...,  N}  . 


•  •  » 


}  is  an  arbitrary  permutation 
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Now  suppose  that  for  each  t  >_  0  there  exist  random  variables 
a^(t)  .^(t) ,  ...»  which  are  H  measurable  and  satisfy  the 

following: 

(i)  ai(t)  >_  0  ,  i  *  1,  N 

(ii)  a^(t)  =0  on  <  t  ,  i  =  1 ,  . . . ,  N 
N 

(iii)  l  ai(t)  *  1  on  T(N)  £  t 

(iv)  The  sample  path  a^C*)  is  measurable. 

The  vector  oi(t)  =  (a^(t) ,  a^(t))  is  the  distribution  program  at 

time  t  ,  that  is,  if  the  load  program  of  the  system  at  t  is  Lg(t)  , 
then  component  i  (i  =  1,2,  N)  sustains  at  t  load  ai(t)Lg(t) 

and  hence  the  failure  rate  of  component  i  at  t  is  61ai(t)Lg(t)  . 
Since  a^(t)  is  a  random  variable,  we  should  interpret  the  failure  rate 
in  this  case  as: 

lim  Ti  P^Ti  €  >  c  +  h]  I  Wt)  =  Siai(t)Lg(t)  w.p.l  . 

h-*0+ 


It  should  be  also  clear  that  the  atoms  of  H (i.e.,  the  elementary 


events)  are  of  the  form 


J 


t^  ,  N(t)  =  kj  and 


since  a^(t)  is  measurable,  it  is  constant  on  these  atoms.  We 

denote  the  numerical  value  of  cx.(t)  on  |t  *t,,...,T  =t 

1  I  "i  1  %  k 

N<0  -  kj  by  at(t  ,  . ,  N(t)  .  k)  . 

This  is  a  generalization  of  the  equal  load  sharing  rule  we  treat  in 


Section  2.1  where  a^(t)  was  there  just 


on  (T  >  t} 


1 

N  -  N(t) 
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All  other  assumptions  about  the  model  remain  the  same. 

We  are  ready  now  to  state  and  prove  lemmas  and  theorems  similar  to 
those  in  Section  2.1. 


Lemma  2.6: 


Proof: 

Clearly  for  A  >  0  : 


P(T(1)  €  (t  ,  t  +  A]  |  Ht)  = 


N 

l  3,a.(t  ,  N(t)  =  0) A  +  o (A) 
jj=l  *  J 

,  as  A  -*■  0  ,  on  {T^  >  t} 


0  ,  otherwise. 


Thus  the  result  follows  by  a  similar  way  as  in  Lemma  2.1. 


Q.E.D. 


Lemma  2.7: 

Let  0  <_  t^  <_  t2  <_  . . .  <_  tfc  <_  t  ,  k<N  ,  and  let  {n^,  ...,  n^}  be 
an  arbitrary  permutation  of  {1,2,  ...,  N}  ,  then  for  x  >  0 
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Proof : 

Similar  to  Lemma  2.2. 

Q.E.D. 

Lemma  2.8: 

{t^}  ,  {n^}  ,  k  ,  as  in  Lemma  2.6  and  let  dt  denote  the  infinitesimal 
interval  (t  ,  t  +  dt]  ,  then  for  x  >  0  , 


Proof : 

Similar  to  Lemma  2.3. 

Q.E.D. 

The  following  generalizes  Theorem  2.1. 

Theorem  2.11: 

Consider  the  simplex  0  <  t,  <  t.  <  ...  <  t..  and  {n,}  is  as  in 

1-2-  -N  j 


Lemma  2.6  then: 
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*(\  £  iil . \  S  i£n)  '  [^s  <'!>%"  •  N<‘l>  •  °>dtl 

'  “  / 

x  II  8  L  (t.)a  It.  ,  T  =  t.  ,  T  =  t .  .  , 

•  o  n.  s  j  a . \  j  n.  1’  n.  .  j-1 

3  “2  j  J  3'  1  3"1 

N(tj)  -  j  -  ljdtj 

()  . 

x  exp  -  <  J  1  ,  N(u)  =  0)Lg(u)du 

(o  j=1 

C2 

+  j  j  ejaj|U  »  Tn  =  C1  '  “  ijLg(U)dU  +  •••  + 

*N  N 

J  £  Vj("  •  \  ■  h . TVl  ■  Vi  •  - N  - 


Proof : 


Similar  to  Theorem  2.1  using  here  Lemmas  2.5,  2.6,  and  2.7. 


Q.E.D. 


As  in  Section  2.1,  we  define  here  the  vector  (W^.W^,  ...,  W^)  ,  where 

T. 

l 

Wi  -  J  ai(u)Lg(u)du  ,  i  =  1,2,  ....  N  . 

0 

Thus,  the  W's  are  the  total  loads  sustained  by  components  1,2,  ....  N  . 
The  following  generalizes  Theorem  2.2: 

Theorem  2.12: 

The  random  variables  W^,  ...,  are  independent  exponentially 
distributed  with  rates  ‘ ‘  respectively. 
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Proof : 


Very  similar  to  the  proof  of  Theorem  2.2. 

Here  is  the  sketch  of  it.  Consider  the  simplex  0  <  t]>  <  c2  1 


<_tN  and  (nr  . . . ,  a  permutation 


of  1,  ...»  N  »  then  on 


L  .  t  T  -  tJ  the  vector  (W . %>  assumes  the 

{\  V  •*  %  *) 


values 


W  *  a  (u  ,  N(u)  =  0)L  (u)du  -  w 
nl  nl  1 


w  *  I  a 

n2  J  n2 


(u  ,  N(u)  »  0)L  (u)du 


C2  v 

J  v(<“>  '  \  *  h  ’  N<U>  '  "I 


l)Lg(u)du  = 


etc. 


Hence  the  joint  density  function  of  W  at  the  point  w  xs  just 


w  .  wn  nl 

n  n«,\  -L 


f_  T  (c_  ,  -  -  •  *  tjj) 

nl  nN 


w  ,  .  .  .  , 


3 (w1>  ....  wn> 


u-1 


9 ( t  ^ ,  •••»  t^) 


Check  that 


£  L.t,J)a»1(tJ  '  \  '  H . V 

C  .  »(tj>  '  J  -  l)  • 
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where  the  first  term  in  the  product  should  be  interpreted  as 

L  (t.)a  (t  ,  N(t..)  =  0)  .  Combine  this  with  Theorem  2.11  to 
s  x  n  ^  X  X 

obtain 


fw(V  • 


N  -3.w. 

n  a  a  3  3 


•  V  “  n  V 

j=l  3 


,  w.  >  0 


The  above  theorem  provides  us  with  most  of  the  results  that  we 
obtained  for  the  equal  load  sharing  rule.  For  completeness  we  state 
them  here  without  a  proof.  Throughout  the  rest  of  this  section  fix 
L  (•)  and  let  (o_(t)  :  t  0}  be  fixed  and  satisfying  conditions 
(2.3).  Note  that  under  these  conditions 


f  N 

I  Lg(u)du  =  l  W  . 
•n  i“l 


Theorem  2.13: 


If 

Ls(*: 

Theorem  2 

.14: 

Fix 

t  > 

.... 

V 

Corollary 

2.9 

Let 

3  » 
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Theorem  2.15 

Suppose 
n  >  1  A  <  3 


Theorem  2.16 

Let  0, 
k 

then 


P<T(N)  ”  c) 


-6  J  Ls(u)du  1f0 


L  (u)du 
s 


1! 


there  exist  0  <  A  <  B  (constants)  such  that  for  any 
<  B  ,  then 


N 

j  L  (u)du  -  l  1/6. 

I  S  1=1  3 


+  N(0,1)  as  N  -►  00 


I  1/8T 

U=1 


1/Bk  ,  k  =  1,2,  N  .  Suppose  that 


-M 


@k  =  0  +  )  as  k  -*■  00  and  0  >  0 


Lg (u) du 


N 


\ 

0/2-  N(0,62)  as  N 


0  . 
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CHAPTER  3 

GENERAL  BREAKDOWN  RULE 

In  this  chapter  we  generalize  the  assumption  about  the  breakdown 
rule.  Whereas,  in  the  previous  chapter  we  consider  only  a  linear- 
power  breakdown  rule,  here  we  allow  it  to  be  quite  general.  We  basically 
adapt  the  model  which  was  treated  by  Coleman  (1957,  1958)  and  more 
recently  by  Phoenix  (1978).  In  all  previous  studies,  the  models 
assumed  parallel  system  of  identical  components.  Here  we  generalize 
it  and  study  systems  of  noniientical  components. 

We  again  consider  an  N-component  parallel  system  which  is  subjected 
to  a  nonnegative  time  dependent  load  program  Ls(')  •  Thus  the  nominal 
load  per  component  is  Lg(t)/N  .  Throughout  this  chapter  at  any  time 
t  >  0  ,  the  surviving  components  are  assumed  to  share  the  load  equally. 

If  we  denote  by  N(t)  the  number  of  components  that  failed  prior  to 
time  t  ,  the  actual  load  carried  by  each  of  the  N  -  N(t)  surviving 
components  at  t  is  Lg(t)/N  -  N(t)  .  Hence  the  actual  component 
load  program  is  a  stochastic  process  which  we  denote  by 

/Lg(t)/N  -  N (t)  on  N (t)  <  N 
L(t)  =  1 

(Ls ( t)  on  N(t)  =  N  . 

Each  component,  therefore,  is  being  subjected  to  L( t)  up  to  its  time 
of  failure  after  which  its  load  is  being  distributed  equally  among  the 
rest  of  the  surviving  components. 

Assumptions  about  the  individual  components:  Component  i  (i  =  1, 
2,  . . . ,  N)  when  it  is  being  subjected  to  a  known  load  program  !(•) 
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has  a  random  failure  time  T  which  is  distributed  according  to 
F^(t  |  !(•))  ,  where 


(3.1) 


F.(t  |  !(•))  =  1  -  e 


01|/  K(l(u))duj 


,  t  >  0 


where  ,  i  =  1,  N  ,  and  K  satisfy: 


(i)  for  each  i  ,  ^  is  increasing,  continuous,  ^(0)  =  0  , 
^i(”)  =  00  • 

(ii)  K  is  positive  increasing,  unbounded  (i.e.,  K(«>)  =  =° 

and  has  continuous  derivative  K'(0  • 


~^i(x) 

If  we  denote  by  G^x)  =  l-  e  ,  i  =  1,  ...,  N  ,  then  Gi 

is  a  distribution  on  the  positive  line  and  F^(t  I  D  =  G±{^S  K0-(u))duj 

If  is  also  differentiable,  then  under  a  known  load  program  1(0  , 

component  i  possesses  a  failure  rate  function: 


(3.2) 


X.(t  |  1) 


4 


K(l(u))du  K(l(t)) 


Note  that  in  this  case  the  failure  rate  function  contains  two  factors 


K(l(u))du| 


'K-IJ  K(l(u))duJ  which  characterizes  the  effect  of  the  load  historv 

and  K(l(t))  which  characterizes  the  effect  of  the  present  load.  As 
in  the  previous  chapter  we  should  emphasize  that  the  components  life¬ 
times  are  stochastically  dependent,  but  only  through  their  failure  rate 
functions  (see  introduction  to  Chapter  2). 
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Our  approach  is  similar  to  the  one  used  by  Phoenix  (1978).  We 
will  show  that  the  system  lifetime,  (under  certain  load  programs)  can 
be  expressed  as  a  linear  combination  of  order  statistics  of  independent 
random  variables  with  known  distribution  functions.  To  obtain  the 
asymptotic  distribution  of  the  system  lifetime  (as  the  number  of  com¬ 
ponents  increase  to  infinity)  we  use  Shorack  (1973),  Stigler  (1974), 
and  some  more  recent  result  of  Wesley  (1977)  which  also  points  out  an 
error  in  the  original  version  of  Stigler.  The  error  in  Stigler' s  paper 
leads  to  an  unjustified  statement  in  Phoenix's  paper  which  we  will  point 
out  and  give  a  correct  version. 


3.1  Preliminaries 


In  this  section  we  review  some  of  the  results  obtained  in  Shorack 

(1973)  and  Stigler  (1974).  Our  interest  in  these  papers  derives  from 

the  fact  that  our  result  leans  heavily  on  them.  Their  results  are 

quite  general  which  is  more  than  we  need,  thus  we  shall  bring  only  a 

modified  version  which  fits  our  case. 

Let  X^  2’  *  * ’ ’  ^n  n  ^e  a  nonnegative  triangular  array  of  row 

independent  random  variables  having  continuous  distribution  fucntions 

n 

F  ,  ,  F  ,  ....  F 

Hj  X  ii)  / 


-  1 

and  F  =  —  J  F  .  and  let  X,  , .  <  X ,  .  < 

n,n  n  n  ^  n,i  (n,l)  -  (n,2)  - 


•••  —  X(n  n)  ke  t^eir  or^er  statistics.  Shorack  and  Stigler  were 
interested  in  the  asymptotic  (as  n  -*■  “)  distribution  of: 


S 

n 


C  .X.  .. 

n,i  (n, l) 


where  C  .  are  known  constants. 
n,i 
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If  we  define  J  on  [0,1]  by 
n 


jn(t) 


C  . 
n,i 


i 


< 


n 


J  (0)  =  Cn  1 
n  n ,  i 


then 


S 

n 


i  n 

1  l 

n  *-• 


J  (i/n)X,  . v 

n  i^1  n'  (n,i) 


Now  suppose  the  following  conditions  hold: 


1.  Let  J  be  continuous  on  (0,1)  and  suppose  J  uniformly 

on  [0  ,  1  -  0]  ,  as  n  -*■  00  for  each  0  >  0  . 

2.  For  some  F  1  ,  F  ^"(t)  ^  F  ^"(t)  for  each  continuity  point  t 

of  F-1  . 

3.  J  ,  Jr  are  uniformly  bounded  on  (0,1)  . 

4.  F  '''(t)  <  D(t)  ,  F  ^(t)  <  D(t)  for  0  <  t  <  1  ,  where 

—  n  — 

D(t)  =  M[t(l  -  t) ]  for  some  5  >  0  . 

Then: 


Theorem  3.1: 

If  conditions  l,  2,  3,  4  hold,  then 

^(S  -  U  )  2  N(0,o2) 

n  n 


_ i 


where 
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1 

Mn  =  I  Jn(u)F'1(u)du 
0 


2 

a 


J(s)J(t)K(t,s)dF_1(t)dF"1(s 


and 


K(s,t) 


sac-  lim 

n 

n-*=° 


Proof : 

Shorack  (1973),  Theorem  3.1. 
a2  can  be  rewritten  as 


a 


2 


00  00 


// 

0  0 


J (F(x))J(F(y))r(x,y)dxdy 


where 


1  n 

T(x.y)  =  lim  -  l  [F  . (x  Ay)  -  F  ,(x)F 
n-*°  n  i=i  n’x  n>1  n> 


Theorem  3.2: 


If  the  above  conditions  hold,  then 

1 

I  J(u)F  ^(u)du  as  n  -*■  <*>  . 
0 


±(y)3  • 

Q.E.D. 
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3.2  The  Result:  (The  Asymptotic  Distribution  of  T^  as  N  -*■  “) 

We  denote  by  T^,T^,  T^  the  lifetimes  of  components  1,2, 

...,  N  respectively  and  by  ^(i)»^(2)>  •••,  T^  their  order  statistics. 
We  assume  that  the  system  is  being  subjected  to  a  nonnegative  load 
program  Lg(*)  and  that  each  component  obeys  the  breakdown  rule  given 

in  3.2.  (That  is,  we  assume  that  ,  i  =  1, ,  N  exists  and  the 

failure  rate  function  under  a  known  load  program  of  an  individual  com¬ 
ponent  is  given  by  3.2.)  To  begin  with,  we  obtain  the  joint  density 
function  of  T^,  . . . ,  T^  under  Lg(*)  using  an  approach  similar  to 
the  one  in  the  previous  chapter.  The  assumptions  of  the  model  described 
before  lead  clearly  to  the  following  facts: 


(FI)  For  each  i  (i  *  1,  ....  N)  , 


PL  <Ti  €  dt  | 
s 


(F2)  If  0-t0<t1<t2<...<tk<t  k  <  j  <_  N  and 

(n^,  ...,  n^}  is  an  arbitrary  permutation  of  {1,  ...»  N}  , 
then 


where 


L (u, •  •  •  9 


Vc) 

N  -  j  ’  Cj  -  C  <  Cj+1  ’ 

J  *  0,1,  ....  k  -  1 

Ls(t) 

-N  -  "k  ’  C  -  Ck 


and  dt^  is  the  infinitesimal  interval  (t  ,  t  +  dt]  . 

The  following  lemmas  are  a  direct  consequence  of  (FI) , 
(F2)  and  their  proofs  are  similar  to  Lemmas  (2.1),  (2.2), 
(2.3)  which  we  omit  here. 


Lemma  3.1: 


VT<» 


Lemma  3.2: 


Let  0  ■  tg  £  t^  <_  t£  •  •  •  2.  ck  <  t  k<N  and  let  n^,  .... 
be  an  arbitrary  permutation  of  1,2,  ...,  N  ,  then  for  x  >  0  we  have 


PL  lT(k+l)  >  C  +  *  *  TnL  "  lV  *'•’  Tnk  =  lk  ’  T(k+1)  >  t) 


N 

'  I 

j-k+1 


n  |  /  K(L(u,t1,  ....  tk))duj  - 

nj(o  Ka(u,ti . ck))du) 


where  L(t,t^, 


•  •  •  * 


t^)  is  defined  in  (F2). 
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Lemma  3 . 3 


{t^}  ,  {ik}  as  in  Lemma  3.2  then  for  x  >  0 


s'  k 


T(k+i)  >  c  +  x  I  Tni  =  cr  —•  \_1  "  Vi  *  T(k)  >  c 


N 

~  I 

i-k+1 


*  ^nj^f  K(^U,tl’  *•*’  Ck-l))d,yK(N  -  k  +  l)dt 

Pn  (  /  K(L(u,t1,  tk))duj  -  ipn  |/  K(L(u,t1, 


•  •  •  *  t<  ) ^ 


We  are  now  ready  to  obtain  the  joint  density  function  of  T^,  . . . , 


Theorem  3.3 


Let  {n^}  ,  { t _^>  be  as  in  Lemma  2,  then 


(T 

A  ni 


.T  E  dt,  i  *  1 

L  \  n,  — i 

s 


. N) 


N 

n 

i=l 


^  KdCu.t^,  ....  tN))d^K(L(ti,t]L,  •••>  t^jd^ 


K(Uu,tr  tN))duJ 


x  e 


where 


(M“> 

IsT^T  ‘ii’'’ 


...»  tM) 


i+1 


i  «  0,1,  ....  N  -  1 


1 . *  ( 


Vu)  u  i  S 


43 


Proof: 

Combine  Lemmas  1,  2,  3  in  a  similar  way  as  in  Theorem  2.1. 

Q  •  E  •  D  • 

Next  we  generalize  Theorem  2.2  to  this  model: 

Consider  ’T(2) »  •••»  ^(n)^  where  L  is  as  defined  above 

and  {T^^}  is  the  order  statistics  and  defined  by: 

Ti 

(3.3)  Qi  “  I  K(L(u’T(l)’  •••»  T(N)>>du 

0 

then  the  following  theorem  is  a  general  version  of  Theorem  2.2: 

Theorem  3.4: 

The  random  variables  Q^,  ...,  QN  are  independent  and  distributed 
according  to  G^,  ...,  respectively,  where 

-'K(x) 

G^(x)  *  1  -  e  x  >_  0  . 

Proof : 

We  prove  it  by  obtaining  the  joint  density  of  Q  ,  . . . ,  Qjj  .  Since 

K(»)  is  (by  assumption)  increasing  and  the  equal  load  sharing  is  applied, 

then  under  L  (•)  >  0  the  mapping  defined  in  (3.3)  is  one  to  one  and 
s 

hence  the  joint  density  of  Q^,  ...,  can  be  expressed  in  terms  of  the 

joint  density  function  of  T^^,  ...,  T^  .  Let  {n^,  ...,  n^}  be  an 

arbitrary  permutation  of  {1,  — ,  N}  ,  if  f  and 

x  »  •  •  • »  x 
nl  "N 
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T  are  the  joint  densities  of  the  Q's  and  T's  re- 

nl  "n 

spectively,  then  f  =f  _  ||j||-1 

9  •••  9  X_  !,•••,!  > 

1  "n  ni  "n 

ii  li  ^  •  •  •  » 

where  |  |J| |  =  3^t  - - y  .  Now  consider  0  <  <  q2  <  ...  <  q^ 

to  find  fQ  q  (qr  •••,  qN)  let  tL  <  t2  <  ...  <  t  be  such 

nl  % 

that. 


1  ■  /  W 


q,  =■  q 


.  •  i  -(K). 


qN  ”  qN-l  +  J 


K(Ls(u))du 


It  is  easy  to  see  that  {t^}  are  uniquely  defined  and  furthermore. 


m 


•m) 


■m 


i 
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Thus: 


V 


-1 


N 

11  ^-i)e 

i*l  i 


1 


<  < 


<2  " 


4N 


Since  this  argument  holds  for  any  arbitrary  {n^}  ,  the  proof  is 
complete. 

Q.E.D. 


Remark: 


The  above  theorem  holds  even  when  the  ij^'s  are  not  differentiable. 
This  theorem,  as  we  pointed  out  before,  is  the  adaptation  of  Theorem  2.2 
to  the  more  general  model  and  its  importance  is  that  it  enables  us  to 
express  T^  —  system  lifetimes  as  a  weighted  sum  of  Q^,  •••* 
the  order  statistics  of  the  Q's  ,  then  use  Stigler's  results  to  derive 
the  asymptotic  distribution  of  T^^  as  N  -*■  ®  .  The  same  approach 
was  used  by  Phoenix  (1978)  to  derive  the  asymptotic  results  for  the 
identical-component  system. 

With  the  function  K(*)  mentioned  in  (3.2)  we  associate  the 
function 


(3.4) 


4>(x,A) 


0<_x<l,A>0 

x  *  1 


9 

From  the  assumptions  about  K(»)  ,  it  follows  that  <Kx,A)  -  $'(x,A) 
is  decreasing  and  continuous  on  [0,1]  . 
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To  illustrate  the  method  we  will  first  consider  the  case  in  which 

the  system  load  program  L  (t)  ■  L  (constant).  Thus  if  T1 ,  . . . ,  T 

s  s  IN 

and  ,  ...»  QN  are  as  defined  in  (3.3)  and  if  T(]_)*T(2)*  •••*  T(^)  ’ 
,  ....  are  their  order  statistics  then  under  Lg(t)  =  NL 

(constant)  we  have: 


T  -  T 

(i+1)  (i) 


T(l)  ' 

♦&L) ' 


Q(i+1)  '  Q(i)  1  =  1*2»  N  ”  1 


thus 


T 


(N) 


Q 


(i)  ’ 


We  define 


(3.5) 


Then 


JN(0)  -  nj*(0,L)  -  <t>(±,L)j  . 


(N) 


H  X  jn(Dq 


i=l 


(i) 


and  furthermore: 

Lemma  3.4: 

uniformly  in  [0,1]  as  N  . 
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i 

k. 


Proof : 

By  assumption,  is  continuous  and  hence  uniformly  continuous 

in  [0,1]  .  By  the  mean  value  theorem  3  t,  £  ( —  ~>^)  such  that 

n|<fr(~  ~  — ,l)  -  ♦(i,L)j  =  — 4> '  ( t ^ ,L)  ,  so  for  each  x  €  [0,1]  3  t  with 
|t  -  x|  <  1/n  such  that  |-4>'(x,L)  -  J^Cx)  |  =  |  — 4> T  (x,  L)  +  <f>'(t,L)| 
which  by  the  uniform  continuity  of  — <}> '  is  uniformly  small  as  n  -*■  00  . 

Q.E.D. 


In  fact 


can  be  rewritten  in  the  following  way 


(N)' 


N  J1  -*'(tN,i’L)Q(i) 


for  some  {tXT  .} 
N,i 


Corollary  3.1; 

{J^}  are  uniformly  bounded  on  [0,1]  . 


Theorem  3.5: 

Suppose  that  the  following  hold: 

(i)  For  some  G  ^  ,  Gn^(t)  G  ^"(t)  for  each  continuity  point 
of  G  ^  ,  where 


G 

n 


l 


i-1 


and  {G^}  are  from  Theorem  3.2. 
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(ii)  G-^  ,  G  ^  are  bounded  above  by  M[t(l  -  t)]  for  some 

n 

6  >  0  ,  M  >  0  . 

Then  under  the  constant  load  program  Lg  =  NL 

^T(T(n)  -  uN)  +  W(0,o2)  as  N  -*•  « 

where 

1 

^  "  /  JN(u)Gn1(u)du 
0 


and 


2 

a 


'  (s,L)4’(t,L)K(s,c)dG~1(c)dG~1(s) 


where 


K(s,t)  =  s  A  t  -  lim  ~  l  Gi(G^1(t))  G1(g"1(s))  . 
n-+«  i=l 


Furthermore  the  asymptotic  mean  of 


is  given  by 


u 


1 

J*  -<p  '  (u,L )G  1(u)du  . 
0 


Proof: 

All  conditions  of  Theorem  3.1  hold. 


Q.E.D. 
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In  his  paper,  Fhoenix  obtained  similar  results  for  the  identical 


components  case.  However,  he  went  on  to  conclude  that  ^(T.  .  =  p) 

In; 

V  2 

-*■  N(0,o  )  .  Unfortunately,  this  is  not  a  direct  consequence  of  the 


above  theorem,  unless  one  shows  that  ET 


(n) 


u  + 


•te) 


as  n  -*■  °° 


Stigler  (1974)  addressed  this  problem  (see  Theorem  4  there)  but  as  it 
turned  out  his  argument  was  not  complete  (see  for  instance  Stigler 
(1979)  or  Wesley  (1977)).  Thus  as  far  as  we  can  tell,  Theorem  2  and 
Corollary  1  in  Phoenix  (1978)  still  need  to  be  justified.  One  possible 
solution  was  provided  by  Shorack  (1972)  (see  Example  1) ,  which  can  be 
stated  as  follows: 

For  the  identical  component  system  (i.e.,  G.  s  G  V  i)  under  a 
constant  load  program  Lg  =  NL  . 


Theorem  3.6: 


(i) 

For  some 

r  >  0 

,  J"xrdG(x)  <  <*>  . 

(ii) 

For  some 

6  >  0 

,  M  >  0  ,  ! 4>  *  (x,L) 

(iii) 

-0- 

exists 

and  continuous  on 

k” 

(x,L)  | 

-  f+: 

£  M[t (1  -  t)] 

— ^5  H - 1-6 

,  r 


then 


^l(T(N)  -  w)  S  N (0 ,o  ;  as  N  +  ® 


where  p  ,  o  are  as  in  Theorem  3.2. 
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Proof : 

See  Example  1,  Shorack  (1972). 


Q . E. D . 


Next  we  show  how  to  derive  similar  results  for  a  more  arbitrary 
load  program.  To  avoid  complicated  expressions,  we  restrict  ourselves 
to  the  power-law  breakdown  rule  (i.e.,  K(t)  =  (Kt)p  ,  p  >  1  and 
suppose  Lg(t)  =  Nl(t)  ,  thus  1 ( t )  is  the  nominal  load  per  component. 
By  the  definitions  of  {T^}  ,  {Q^^}  we  have 


Q 


(1) 


[K1 (u) ]Pdu 


Q(i+1)  “  Q(i) 


1. 


N  -  1 


thus 


(l(u))Pdu 


K  PQ 


(1) 


T(i+1)  . 

j  a(u))Pdu  -  -  w 

T 


ana 
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Define 


ine:  <f(x,k)  =  ( 1~;  — )  . 


Theorem  3.7: 


If  Che  conditions  of  Theorem  3.2  hold,  then 


(N) 

J 


(l(u))  du  -  p 


V  2 

-*■  W(0,o  )  as  N  •>  00 


where  u„  ,  a  are  as  defined  in  Theorem  3.2.  Suppose  1(0  >  0  , 
N 

and  if 


-/ 


h(t)  =  J  [  1  (u)  ]  du  ,  t  _>  0 
0 


then  h  1  is  well  defined  and 


Corollary  3.2: 


as  N  00 


where  u  ,  ,  a  as  in  Theorem  3.2. 


Proof : 


The  above  is  a  special  case  of  Slutsky  theorem. 


Q.E.D. 
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